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Using heavy-baryon chiral perturbation theory, we have calculated all the diagrams 
up to two-loop order which contribute to the <S-wave pion self-energy in symmetric nuclear 
i -^h ' matter. Some subtleties related to the definition of pion fields are discussed. The in-medium 

j—h ' pion mass turns out to increase by only (5 ~ 6) per cents at normal nuclear matter density, 

^ ' and this result is model independent. 
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The role of pionic interactions in nuclear medium has long been an important issue in 
nuclear physics. Recently, the interest in the S-wave pionic interaction in nuclear medium 
has reemerged due to the experimental observation of deeply bound pionic states in Pb 
at GSI |, §. It is to be noted that there is a wide variance among the theoretically 
deduced values of the effective pion mass m* ; for symmetric nuclear matter, m* is found 
to increase by 7 % in Ref . |2| , but by 20 % in Ref . || . These results are obtained relying on 
the phenomenology related to the <S-wave pion-nucleus optical potential. It is well-known 
that threshold isoscalar S-wave pion-nucleon interaction in nuclear medium, though very 
weak, is more repulsive than in vacuum |4|. Noting that the leading Weinberg- Tomozawa 
term vanishes and the pion-nucleon sigma term is partly canceled by the effective range 
term, higher order corrections to the one-body contribution is substantial ||, and also the 
two-body contribution is expected to play a significant role; this latter aspect has been 
considered in explaining the threshold S-wave pion-nucleus interaction [||, ^, |7j]. 

There are approaches for the effective pion masses based on chiral symmetry J8|. 
In this letter we reexamine the problem in symmetric nuclear matter in the framework of 
heavy baryon chiral perturbation theory. All the contributions up to 0{Q I ) and the one- 
body 0(Q 8 ) contribution in the pion self-energy will be calculated. In addition, we will 
discuss some subtleties of the problem related to the definition of the pion fields. 

The relevant quantity for the effective pion mass m* in nuclear matter we are in- 
terested in is the in-medium self-energy of pion, II(u; 2 ), which is defined with the pion 
momentum = (uj, 0). To evaluate the self-energy systematically, we employ Weinberg's 
formalism Q of heavy baryon chiral perturbation theory (HBChPT) , which has been quite 
successful in understanding threshold pion-nuclear interactions [||, [l0|, |ll , |12 ] . In Weinberg's 
scheme, Feynman diagrams are characterized by the order Q u , where Q is the typical size of 
the momenta involved and/or pion mass, which is regarded as small compared to the chiral 
scale A ~ 1 GeV. In this scheme, the Fermi momentum kp is thus counted as of order Q. 
For diagrams for the pion self-energy, the chiral index v is given as 

v = 2 + 2L + ^2 v i, v i = d i + y - 2 (!) 

i 

where L is the number of the loop, di is the the number of derivatives and/or pion masses, 
and rii is the number of nucleon fields attached to the i-th vertex. In HBChPT, the La- 
grangian consists of pions and nucleons, while all other degrees of freedom are integrated 
out, 

£ eff = £ + £ 1 + £ 2 + ---, (2) 

where the subscripts are the chiral index V{. Terms relevant to computing the S-wave pion 
self-energy are 

£ = B[ivD + 2ig A S ■ A] B - ^ C A {BT A B) 2 + /*Tr (WiA M + *±^j , (3) 
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The low-energy constants ci, C2 and C3 are to be fixed from experimental data; at present 
there are several sets of values for the parameters determined by fitting low-energy pion- 
nucleon data||: 

(ci,c 2 ,c 3 ) = (-1.27 ±0.12, 3.23 ±0.19, -5.93 ±0.08) GeV" 1 (Fitl), 
(ci,c 2 ,c 3 ) = (-1.47 ±0.09, 3.21 ±0.11, -6.00 ±0.03) GeV" 1 (Fit2), 
(ci,c 2 ,c 3 ) = (-1.53 ±0.18, 3.22 ±0.25, -6.19 ±0.09) GeV" 1 (Fit3). (6) 

Figure 1 shows the nuclear diagrams we compute for the pion self-energy. Note that 
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Figure 1: Feynman diagrams that contribute to the pion self-energy. Solid lines are nucleons, and 
dashed lines are pions. The rectangle in (a) denotes the whole set of diagrams (including loop 
corrections) that contribute to the ttN scattering amplitude. 



nucleon loops in free space are forbidden in HBChPT, and all the nucleon lines (in solid line) 
drawn in Fig. 1 run only up to k p, the Fermi momentum. Throughout this work, we employ 
the Fermi gas approximation for the wavefunction of nuclear matter. We acknowledge that 
the diagrams in Fig. 1 have been studied in chiral models since long time ago, see, for 
example, Ref. [|l3| . 

The leading order (LO) contribution to the 5-wave pion self-energy is from the one- 
loop graph with the Weinberg- Tomozawa term coming from Cq, which vanishes for sym- 
metric nuclear matter. If such a term survives, it would be 0(Q 4 ). The next-to-leading 
order (NLO) contribution comes from the one-loop graph with the ttttNN vertices from £\, 
which is nothing but the pion-nucleon sigma term. Before going further, we note that all 
the diagrams with a one-nucleon-loop only, called one-body contributions, are related to 
the pion-nucleon scattering amplitude. Particularly, all the contributions at u = m^ can 
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be replaced by the isoscalar irN scattering length, a n N 
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^ 2L . Note that the above equation is correct to all order, as far as the one-body 
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contributions are concerned. The scattering length has been calculated within HBChPT 
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In the r.h.s. of the above equation, the first term is from 
NLO, the second from the next-to-next-to-leading order (NNLO), and the ellipsis denotes 
higher-order contributions. At NNLO, we also have genuine two-body contributions which 
are drawn in Fig. l(b — d). Including all of them, the pion self-energy up to NNLO or up 
to O (Q 6 ) reads 
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and C = | f° r the pion field defined by eq.(|5|). Here, £ is an off-shell parameter (that is 
effective only when the pion is off-shell), which depends on the definition of the pion field. 
For example, the canonical pion field introduced by Weinberg [||, 14 1 ^ corresponds to £ = 3. 



In fact, the value of £ is not fixed by chiral symmetry and can be any arbitrary number 
[14]. As another interesting definition of the pion field, we discuss here the background field 
method (BFM) [ffq]. In BFM, the covariant quantities are defined as 



£r = exp 
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where 7r c i is the pion field that appears as external legs, while the ttq is for the internal pion 
lines. Often the BFM removes some of artificial contributions and simplifies the calculation. 
In our case, the contribution of Fig. (lc) vanishes with BFM; with the canonical pion field 
of eq.(|5|), Fig. (lc) is not zero but canceled out by a part of the contribution of Fig. (Id). 
For physical observables where all the external particles are on mass shell, however, BFM 
gives exactly the same results as any other definition of pion field does. The value of £ 
obtained by using BFM is £ = 4. 
#1 The covariant derivative of the pion field we are referring to has the form <9 M 7f/(l + 7r 2 /(4/^)). 
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Figure 2: The ratio m*/m„- as a function of p/po- The upper shade is obtained with the canonical 
[/-field (C = |), while the lower shade with the BFM (£ = 4). In each shade, the upper (lower) 
bound corresponds to Fit 1 (Fit 3). Fit 2 gives almost the same results as Fit 3. 



Now we are in a position to discuss the in-medium pion mass m*. The conven- 
tional way of determining the in-medium mass is to find the zero of the inverse of the pion 
propagator, w 2 — m 2 — U(w 2 ), or 

mf=ml + U(mf). (12) 

But since II(m* 2 ) with m* ^ m, depends on £, the resulting m* also depends on the 
off-shell parameter £. The results are depicted in Fig. 2 with two popular choices of £; 
the upper shade is for the case with £ = 4/3 (the canonical field of eq.(H)), and the lower 
shade is for £ = 4 (BFM). The results with £ = 3 lie in between the two. In each case, the 
parameter set Fit 1 gives the biggest m*, and Fit 3 results in the lowest value. At normal 
nuclear matter density, the ratio m^/m^ is 1.20 — 1.32 with Q = 4/3, and 1.12 — 1.15 with 
£ = 4. Thus, apart from the non-trivial dependence on the parameter set of low-energy 
constants, there is quite substantial ^-dependence. We note that in the chiral limit = 0, 
however, m* remains zero in all the cases. 

An interesting observation is that, when the nuclear density reaches around p = 
(1.3 ~ 1.7) po with C = (4/3 ~ 4), where po is the normal nuclear matter density, the in- 
medium pion mass decreases abruptly below the pion mass in the vacuum. This is caused 
by the fact that when the density increases, the value of the IT(u; 2 ) = -^^H(lo 2 ) is close to 
one and makes the solution of eq.(|T2|) unstable at the density mentioned above. 



The off-shell ambiguity may indicate that the effective pion mass defined by eq.(12) 
is not a physical observable, which should be independent of the definition of the pion field. 
Indeed the effective pion mass is measured only indirectly in experiments with the help of 
optical potentials. A thorough study on this subject, which should be quite important, has 
not been done in this work and will appear in a separate work. 

In this letter, however, we introduce a new way to extract the effective pion mass 
unambiguously. So far the pion self energy has been obtained to a certain order u, and 
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then treated non-perturbatively in obtaining the effective mass through eq.(12). To avoid 
the above mentioned ambiguities, we propose to get the effective mass itself up to order v. 
The new procedure gives us exactly the same results as the previous method up to u, all 
the differences residing in higher orders. Particularly if the gap (m* 2 — m 2 ) is of order of 
Q 2 as the naive counting indicates, there is absolutely no difference. In actual cases, the 
gap is much smaller than in the naive counting, 
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0(Q 5 ), (13) 



since the non- vanishing-leading-order (or NLO) of the pion self-energy is O (Q 5 )- Inserting 
this into eq.(|]) with to 2 = m* 2 , we see that all the (-dependent terms in II(m* 2 ) are of 
0(Q 9 ), which is beyond the scope of this work and will be neglected. In other words, there 
is no off-shell ambiguity in IT(m* 2 ) up to 0(Q S ). Furthermore, II(m* 2 ) is the same as 
n(m 2 ) up to O (Q 7 ), that is, n(r< 2 ) = n(m 2 ) + O (Q s ). Expanding n(r< 2 ) based on our 
"detailed" power counting rule, 

oo 

nK 2 ) = ^nMK 2 ), (14) 

our two-loop order calculation determines the first two leading order contributions unam- 
biguously: 

U^(m* n 2 ) = ^(I + mKat-^, (15) 
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We remark that Il( 5 )(m* 2 ) in eq.(|l5D includes all the one-body contributions at uj 2 = m 2 
up to infinite order, since the experimental scattering length a^N contains all the higher 
order contributions as well as the leading O (Q 2 ) contribution. We can also show that there 
is no 0(Q ) contribution, 

n( 7 V; 2 ) = o, (17) 

by the following arguments. All the one-body contributions are already absorbed into II^ 5 ) 
at uj 2 = m 2 , and the difference II(m* 2 ) — H(m 2 ) is O (Q 8 ), which can be read from eq.@. 
The two-body contributions at O (Q ) are two-loop diagrams which contain one vertex from 
jCi; however, they are all isovector and therefore cannot contribute to the pion self-energy in 
symmetric nuclear matter. Contributions involving three or more nucleons begin to appear 
only at 0{Q 8 ). 

At O (Q s ) there are one-, two- and three-body contributions, 

n( 8 V; 2 ) = <2 body (m; 2 ) + n^ body « 2 ) + n® body (< 2 ). (is) 
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The one-body contribution at this order comes from the w-dependence in n(u; 2 ), 

<2 body « 2 ) = Mi + mKat ^? (0) ■ (19) 

The two- and three-body contributions at this order are not discussed in this work. 
1.1 -i r 
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Figure 3: The cumulative graphs for m^/m^ vs. p/po- The dotted, dotted-dashed and solid lines 
stand for the results up to O (Q 5 ) , O (Q 6 ) and O (Q 8 ) . The result up to O (Q 8 ) does not contain 
the two- and three-body contributions at this order, see the explanation in the text. 

The numerical results are given in Fig. 3. The O (Q 5 ) (dotted line) contribution is 
tiny because of the smallness of the ir — N scattering length. The two-body contribution 
embodied in O (Q & ) (dotted-dashed line) gives the main contribution, about 6 % increase 
of the in-medium pion mass at normal nuclear matter density. As mentioned above, there 
is no O (Q 7 ) contribution. Just for discussion, we have drawn O (Q 8 ) results (solid line) 
without the two- and three-body contributions. This (partial) O (Q s ) contribution at this 
order is not significant, which is again due to the smallness of the ir — N scattering length. 
For the same reason, the dependence on the parameter set up to this order is tiny and 
invisible in the figure. 

By expanding II(m* 2 ) with the detailed power counting rule which considers the gap 
(m* 2 — m 2 ) to be of order Q 5 , we find that the in-medium pion mass increases only very 
mildly as the nuclear matter density increases. At normal nuclear density, m* is larger than 
only by (5 ~ 6) %. This result has little dependence on the parameter set used, and is 
independent of the off-shell parameter Q. Our results support the phenomenological study 
given in Ref.Q. The complete calculation up to O (Q 8 ) is feasible and involves computing 
three-body contributions as well as one- and two-body contributions, and it would play 
more important roles in higher density regions. 

Let us make concluding remarks. We have considered two schemes for computing the 
effective pion mass in symmetric nuclear matter based on chiral perturbation theory: a navie 
power counting scheme (eqs.(|7Hl~0|, |l2| ), Fig. 2) and a more detailed power counting scheme 
(eqs. (15-19), Fig. 3). The discrepancies in m* 2 are of order of 0(Q 9 ) or higher order. We 
believe the latter is more meaningful than the former, because the latter is independent of 
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the offshell parameter £. While of higher order than considered, the possible dependence on 
the off-shell parameter at 0(Q 9 ) remains to be addressed. Especially the radical differences 
between the two schemes (Fig. 2 and Fig. 3) are quite alarming. Thus our results up to 
NNLO should be taken with caution. To be more definite, a calculation up-to-0(Q 9 ) would 
be required. 

In overcoming the problem caused by the off-shell ambiguity, the idea of Brown- Rho 



(BR) scaling [17] may shed some light on this problem. If we follow the BR scaling, the the 
parameters of the pion self-energy should be replaced by the in-medium effective values. 
Let U*(w 2 ) be the BR scaled pion self-energy. Then the equation for the m* reads 



+ rr 



* 2\ 



with 



lT(m; 2 ) 



77U 



, * 2 



2k 3 

-4vr(l + ukl/v— | + 



*1 



2kt 



9a 



J F {ml\ 



(20) 



(21) 



We note that there is no (-dependence, and these equations are essentially the same as 



eqs.(15 1S) except that all the parameters are replaced by their in- medium counterparts. 
This is because in nuclear medium the effective on-shell condition is u> = to* and II*(m* 2 ) 
is effectively an on-shell quantity. However, in order to perform a numerical analysis, we 
need the value of a* nN as well as other effective parameters, which requires more studies. 

On finalizing this work, we have learned of the existence of a similar work by Kaiser 
and Weise(l^]. Their work is more extensive than ours in that it includes pions and kaons 
in isospin asymmetric nuclear matter but the effective pion mass is simply assumed to be 



ml + n(m 2 ). 
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